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Single point spacecraft observations of the turbulent solar wind flow exhibit a characteristic non- 
axisymmetric anisotropy that depends sensitively on the perpendicular power spectral exponent. 
For the first time we use this nonaxisymmetric anisotropy as a function of wave vector direction to 
test models of magnetohydrodynamic (MHD) turbulence. Using Ulysses magnetic field observations 
in the fast, quiet polar solar wind we find that the Goldreich-Sridhar model of MHD turbulence is 
not consistent with the observed anisotropy, whereas the observations are well reproduced by the 
"Slab + 2D" model. The Golderich-Sridhar model alone can not account for the observations unless 
an additional component is also present. 
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In-situ satellite observations of the solar wind show a 
broad spatio-temporal range of plasma fluctuations [l| 
and a high effective magnetic Reynolds number 0, [J. 
This suggests a nonlinear cascade consistent with turbu- 
lence 0, Q . In-situ observations have hence been used 
extensively to test the predictions of Magnetohydrody- 
namic (MHD) turbulence In the solar wind the 
background magnetic field imposes a preferred direction 
in the plasma resulting in an inherent anisotropy in the 
turbulent fluctuations 12 , [l3| ■ This anisotropy is fun- 
damental to the understanding of cosmic ray scattering 
Q and the non-adiabatic rate of cooling observed for the 



solar wind [14 1 



A pioneering anisotropic model of MHD turbulence is 
that of Goldreich and Sridhar [l^ (GS model hereafter). 
In the GS model, the turbulence is dominated by the 
cascade perpendicular to the local mean magnetic field 
which follows Kolomogorov phenomenology with a one- 
dimensional perpendicular energy spectrum predicted to 
be E{k±) ~ k^^^^. A "critical balance" between the lin- 
ear and nonlinear time-scales results in a one-dimensional 
parallel spectrum of ii^(fc||) ~ for the GS model. Here 
k± and fcy are the components of wave vector k perpen- 
dicular and parallel to the local magnetic field direction. 
The energy spectrum is intrinsically axisymmetric w.r.t. 
the local magnetic field. There is significant numeri- 
cal and observational support for the predictions made 
by the GS model. This support includes high resolu- 
tion direct numerical simulations (DNS) JJJ and recent 
solar wind observations that confirm predictions of the 
anisotropy of the spectral exponents [1, [13, 18 1. Para- 



doxically, the GS model is inefficient in scattering cosmic 
rays due to the inherent strong anisotropy {k± ^ fc||) 
at small scales, which cannot account for the observed 
scattering characteristics 0- 

A possible resolution of this inconsistency is that there 
is an additional 'slab' component of solar wind fiuctua- 
tions that is independent of the perpendicular cascade. 



consisting of fluctuations with wave vectors parallel to 
the local magnetic field. There is considerable observa- 
tional evidence that in the solar wind turbulence there is 
a combination of perpendicular and parallel wave num- 
bers. Classically this is seen in the magnetic field cor- 
relation function which, assuming axisymmetry, gives a 
'Maltese cross' pattern [^, [l^. This result was inter- 
preted as a superposition of two-dimensional (fluctua- 
tions with k purely perpendicular to the local magnetic 
held) and slab components [3] (S2D model hereafter). 
The observed anisotropy of the spectral exponents that 
support the GS model in the solar wind do not exclude 
the S2D model [I,[l3,[l3|- is thus an open question as 
to whether the GS model operates as the sole turbulence 
process in the solar wind. 

Belcher and Davis used Mariner 5 observations to 
investigate the anisotropy of solar wind magnetic fluc- 
tuations in the low frequency (energy containing) and 
inertial intervals. They observed that the power in the 
fluctuations is ordered both w.r.t. the average magnetic 
fleld and the solar wind flow direction. If the solar wind 
turbulence is described by the S2D model with the same 
spectral exponents for both components, this observed 
nonaxisymmetric anisotropy can be explained as a sam- 
pling effect [10] related to the Taylor hypothesis [2^. Im- 
portantly, if the parallel and perpendicular components 
have distinct power spectral exponents, as in models such 
as GS, then the Taylor hypothesis can quite generally 
relate the exponents to observed nonaxisymmetry 21 1. 
This suggests a new quantitative test for theoretical pre- 
dictions, since the observed nonaxisymmetry does not de- 
pend solely upon the power spectral exponents of E(k±) 
or E{k\\), but as we shall see, is sensitively dependent 
upon the transition from k± to fcy described by the mod- 
els. In this Letter we will use this idea to test the pre- 
dictions of the GS model and the S2D model against so- 
lar wind observations. We find that the S2D model can 
fit the data whereas the GS model cannot. In addition 
to critically balanced turbulence some other component 
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must be present in the solar wind fluctuations to account 
for the observations. 

Theories of turbulence predict how the energy in the 
fluctuations vary with scale. For fluctuations in a vec- 
tor field, the Fourier transform of the two point corre- 
lation matrix Rij{r) = {SBi{yL)6Bi{x + r)) then defines 
a spectral energy density tensor Pij{k) [ll|, [21 1 which 
captures the full anisotropy of the fluctuations. Single 
point observations in the flow cannot isolate k uniquely 
and instead give a reduced one-dimensional spectral ten- 
sor Pij{f). From Taylor's hypothesis [i^], the measured 
one-dimensional tensor Pij{f) is the integral of Pij(k) 
over the plane k • V^^ = 27r/ in k space, where Vs^u is 
the solar wind velocity, i.e. we observe: 



(fkPj{k)6(27rf - k-Vs 



(1) 



The observed spectral tensor Pij{f,9) is in general 
nonaxisymmetric, with dependence on both / and the 
angle 6 between the local magnetic field and the solar 
wind flow velocity. This is the case even if the un- 
derlying turbulence phenomenology is axially symmetric 
w.r.t. the local magnetic field direction, as is the case 
for both the GS and S2D models. A natural co-ordinate 
system [l^, U, 21 [ is to project the magnetic fiuctua- 
tions onto a local, scale dependent mean field direction 
e^{t, f) = B{t, /)/ |B(i, /)| where B(i, /) is the scale de- 
pendent local average field. If the bulk fiow velocity unit 
vector is V then the other two perpendicular unit vectors 
of the set are: 

e X V 

exitj)^— Tj, ey(i,/) = X e^,. (2) 



e. X V 



We will focus on the power ratio in this coordinate sys- 
tem: 



R(0 f) — -^-^ 
PyyiOjy 



(3) 



We now calculate ([3]) for the GS and S2D models and 
compare with observations. 

GS model: We will assume that for the GS model all 
parallel components of the spectral tensor vanish so that 
Pzi = 0, consistent with other studies [l^l- The perpen- 
dicular components for the GS model are related to the 
power tensor in the following manner: 



pGS 



(k) 



(4) 



where Hij — Sij — k±ik±j/k'j^ and from Ref. 15 1: 



E{k) - Ck 



£2/3^1/3 /fc|,Ll/3' 



.10/3 



1.2/3 



(5) 



Here Ck is the Kolmogorov constant, L is the character- 
istic injection scale, e is the energy dissipation rate and 
5 is a positive symmetric function related to the scaling 
between k\^ and k±, where g{0) = 1 and g{z)dz — 1 
15|. The theoretical prediction of the GS model is 
determined by substitution of Eqs. ([1]), Q and ([5]) 
into Eq. ([3]) and numerical integration. We find that 
the result is not sensitive to the functional form of the 
scaling function g and we use the exponential function 
[5 = exp(-Li/3|fc|||/A:f )]. 

The GS model, as discussed here, assumes balanced 
MHD turbulence with equal power in the Alfven waves 
travelling in both directions along the local magnetic 
field. Turbulence in the fast solar wind is known to be 
imbalanced with more power in the Alfven waves propa- 
gating outward from the sun than toward it [1]. However, 
an imbalanced extension predicts the same spectral 
exponents of both dominant and subdominant waves as 
well as the same scaling relation between fcy and k± as the 
balanced GS model. Thus, given isotropy of the turbu- 
lence at the injection scale L, imbalance of the turbulence 
does not affect the ratio ([3]). 

S2D model: We follow the prescription of Ref. [lo| . 
except that we specify distinct spectral exponents qs and 
q2D for the slab and 2D components as is observed. In 
this Letter, we will consider variation in q2D, as R{9,f) 
is sensitive to the perpendicular spectral exponent. For 
sirnplicity a constant value of qs = 2 is used throughout 
[l,|2lj. In this case the ratio i?s2(^', /) predicted by the 
S2D model is 



A ( fL y I ( fL y 

V Vs„ cos 9 y ' \^ V,^, sin B ) 



where L is the injection scale (we assume that both com- 
ponents have the same injection scale) and Ac is a con- 
stant that characterizes the energy ratio of slab and 2D 
components. Since qs ^ q2D the relative power of the 
components i s / dependent, as is the ratio Rs2{6,f)- 
Observations [lOl give a 1 : 4 energy ratio between slab 
and 2D components and we will fix this at the injection 
scale L, which fixes the constant Ac « 0.5. 

We plot these theoretical predictions in Figure 1, for 
9 — [5°, 85°], as co-ordinate system ([2]) becomes unde- 
fined as 6* — > 0° . The ratio R is plotted using normalized 
frequency F — 2tt f L /Vsw and we show two different cuts 
through the surface R{9, F); for F = 36 (upper plot) and 
F — 3 (lower plot) in the left panels and for 9 = 25° (up- 
per plot) and 6* = 45° (lower plot) in the right panels. 
The S2D model is indicated by red symbols on the plot, 
with circles for Kolmogorov q2D =5/3 and diamonds for 
Iroshnikov-Kraichnan q2D = 3/2. The GS model predic- 
tion is indicated by the black squares. 

For underlying turbulence that is axisymmetric, with 
E{k±_) - k~''^ the observed R{9,f) 7_l for 9 90° 
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Figure 1: Left panels - 6 dependence of the ratio R{6, F) at two normalised frequencies: _F" = 36 (upper plot) and F — 3 (lower 
plot). Right panels - frequency dependence of the ratio R(9,F) aX 6 = 25° (upper plot) and 9 — 45° (lower plot). Observations 
are shown by blue stars, S2D model is in red with circles for Kolmogorov and diamonds for Iroshnikov-Kraichnan perpendicular 
scaling. A solid red line indicates the S2D model with the perpendicular scaling exponent from the data of q2D ~ 1.59. All 
S2D models shown here use qs = 2. The GS model is shown by black rectangles. Middle panels - Number of samples to form 
7?, where white = 1800 and black = 36 - this varies significantly with both frequency and angle. 



for all of these curves [21| . It can be seen in the left pan- 
els of Figure [1] that the predicted R{9, f) for the GS and 
S2D models are distinct for intermediate values of 9 in 
the transition R{e = 0° ^ 90°, /). The GS model form is 
concave, whereas the S2D model is strongly convex. This 
provides a strong test against observations provided the 
statistical variability is smaller than the difference be- 
tween the predicted curves. Importantly, the model pre- 
dictions are maximally distinct for intermediate angles 
(i.e. 9 ^ 20° — 40°) and as we shall see, this is where the 
observations tend to be more statistically significant as 
there more samples. 

Observations: We use magnetic field data at 1 second 
resolution for an interval [day 91-146, 1995] of fast solar 
wind observed by Ulysses. In this interval Ulysses moved 
from a heliographic latitude of 21° to 58° and a radial 



distance of 1.36 to 1.58 AU. This long quiet interval of 55 
days is needed to obtain good statistical coverage across 
9. This interval is of fast solar wind with an average 
fiow speed of 756 km/s and average plasma parameters: 
magnetic field |B| ~ 2.9 nT, ion plasma /3 ~ 1.35, ion 
plasma density ~ 1.23 cm~^, ion temperature Ti ~ 21 
eV and Alfven speed ~ 56 km/s. 

We use the continuous wavelet transform (CWT) with 
a Morlet wavelet to resolve vector fluctuations (5B(i, /) 
in time t and frequency /. The scale dependent local 
field B(t, /) is calculated via the convolution of a Gaus- 
sian window as outlined in Ref. [2l[ to obtain the unit 
vector direction ez{t,f ) = B(t, /)/ |B(t, /)|. In practice, 
the polar wind seen by Ulysses is stable and within 3° of 
the radial direction over the entire interval under study, 
so we replace V by the radial unit vector here, which 
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Figure 2: The surface /) is siiown in colour. The black 
solid lines indicate the cuts of /) shown in Figure [1] The 
dashed line shows the start of the inertial range, F = 1. The 
dash-dot lines show contours of subinterval sample size. The 
colour bar shows the ratio R{Q, f) — Pxx{6, ,f)/Pxx{d, f) 



corresponds exactly with the coordinate system used by 
Ref. 12]. The magnetic fluctuations are then projected 
to basis ([U to determine the power for different compo- 
nents of Pii(t,f) as a function of frequency /, such that 
Pii{t, /) = 2ASBi{t, /)^, where A is sampling time of the 
data. 

For comparison between the models and the observa- 
tions the data is normalized to = 1 at the beginning 
of the inertial range. This corresponds to a time-scale of 
25 minutes. Ratio and a measure of the statistical 
variability for four specific cases of d and F are indicated 
in Figure [T]by the blue stars. The method for estimating 
the statistical variability in these observations is given 
below. Given this variability, we can immediately see 
that the observations correspond quite well to the S2D 
model whereas the GS model fails to predict the obser- 
vations. Thus, the GS model cannot reproduce the ob- 
served inertial range fluctuations in the fast solar wind 
in the absence of some additional, scaling component of 
fluctuations. For completeness, we use the data to obtain 
q2D — 1-59, this gives an S2D model result shown by the 
solid red line. 

The full surface R{6, f) for our interval is shown in 
Figure [5] This shows how the nonaxisymmetry between 
the perpendicular directions depends on the sampling do- 
main. The surface is constructed by subdividing the en- 
tire interval into 1 day subintervals. Pxx{t, f)/ Pyy{t, f) 
is binned according to 6 for each of these subinterval. 
As the distribution in each bin has a form close to log- 
normal the best measure of the subinterval average, R, is 
the geometric mean [l^. In Figure [5] the contours indi- 
cate the number of samples in each R, thus the contours 



may be interpreted as confidence contours. Each bin of 
the surface R{9, f) has 55 realisations of R. These re- 
alisations are used to calculate the statistical variability 
of the distribution R{0, /) by calculating the median and 
interquartile range of the 55 values of R. Figure [5] shows 
the median value for each bin and the stars in Figure 
[1] show the median value with statistical variability esti- 
mated to 99% certainty indicated by the errorbars. Ob- 
servationally, since cos0 = ez-ea — B(t, /)/ |B(i, f )\-eR 
any uncertainty in B(i, /) is greatest as — 0°. We es- 
timate that an uncertainty of 1% on B(i, /) translates to 
S9 - 8° in the 0^0°- 10° interval. Thus the observed 
values at small 9 are unreliable. 

Much of the observational support for GS has re- 
lied upon from the determination of spectral exponents. 
Whilst this is not an unreasonable test for the theories 
and models of plasma turbulence it is not complete and 
may lack uniqueness. Our work highlights the need for 
other distinct methods to test these predictions. Such 
tests need to probe the full three-dimensional energy 
spectrum as the kinematic power ratio used in this Letter 
begins to do. 
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